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Abstract. The CDF experiment is successfully collecting data from pp collisions at the Tevatron in 
Run II. As the data samples are getting larger, systematic uncertainties due to the measurement of the 
jet energy scale assessed using the calorimeter simulation have become increasingly important. In 
many years of operation, the collaboration has gained experience with GFLASH, a fast parametriza- 
tion of electromagnetic and hadronic showers used for the calorimeter simulation. We present the 
performance of the calorimeter simulation and report on recent improvements based on a refined in 
situ tuning technique. The central calorimeter response is reproduced with a precision of 1-2%. 
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1. INTRODUCTION 

Since the start of Run II in 2001, the Collider Detector at Fermilab (CDF) [1] has col- 
lected over 1 fb _1 of data of pp collisions at 1.96 TeV center-of-mass energy. In a variety 
of resulting publications, the simulation of the calorimeter has proved to be a crucial el- 
ement for the precision measurement of physical observables, like the mass of the top 
quark, since it appears as one of the keys to control the jet energy scale systematics. 
The CDF calorimeter simulation is based on GFLASH |2j], a FORTRAN package used 
for the fast simulation of electromagnetic and hadronic showers. It is embedded in a 
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GEANT3 framework ||3|] as part of the whole detector simulation and has various advan- 
tages w.r.t. the detailed GEANT shower simulation: In CDF it is up to 100 times faster, 
and it can be flexibly tuned. The calorimeter simulation was initially tuned to test beam 
data and has been improved due to a steadily refined in situ tuning using samples of 
isolated charged particles, which recently became available over a remarkably extended 
momentum range up to 40GeV. Here we report on the GFLASH performance effective 
for current CDF physics publications and on ongoing improvements contributing to the 
ambitious Run II physics program. 

2. CDF CALORIMETRY 



CDF is a general-purpose charged and neutral particle detector with a calorimeter and 
a tracking system. The calorimeter has a central and a forward section subdivided into 
a total of five compartments (Fig.[T]): the central electromagnetic, CEM Mm, the central 
hadronic, CHA Q, the plug electromagnetic, PEM [6], the plug hadronic, PHA |@], 
and the wall hadronic, WHA [5]. The calorimeter is of sampling type, with lead (iron) 
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FIGURE 1. Le/f." Quadrant view of the CDF calorimeter showing the electromagnetic (CEM, PEM) 
and hadronic compartments (CHA, WHA, PHA). Right: Sampling structure and energy resolutions for 
the central (CEM, CHA) and the plug part (PEM, PHA). 



absorbers for the electromagnetic (hadronic) compartments, scintillating tiles and wave- 
length shifters. It is subdivided into 44 projective tower groups, each group made of 
24 wedges (partially 48 in the plug part) circularly arranged around the Tevatron beam 
and pointing to the nominal event vertex. The central and plug part together cover the 
pseudorapidity range |rj| < 3.6. The calorimeter encloses a tracking system consisting 
of a vertex detector and a cylindrical drift chamber, both situated within a 1 .4 T solenoid 
magnet. It provides a precise measurement of single charged particle momenta serving 
as an energy reference for the calorimeter simulation tuning. 



3. GFLASH IN A NUTSHELL 



The CDF simulation uses GEANT to propagate particles from the main interaction point 
through the detector volume. A shower in GFLASH is initiated when a particle undergoes 
the first inelastic interaction in the calorimeter. GFLASH treats the calorimeter as one 
effective medium using GEANT geometry and material information. GFLASH is ideal for 
calorimeter modules which may have a complicated but repetitive sampling structure. 
Given an incident particle energy, E mc , the visible energy in the active medium, 



d£ vis (r) =E mc m 
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is calculated according to the sampling fractions for electrons (e) and hadrons (h) 
relative to the sampling fraction for minimum ionizing particles (m), taking their relative 
fractions c em and Chad of energy deposited in the active medium into account. / e m( r ) 
and /had( r ) are electromagnetic and hadronic spatial energy distributions of the form 



/(r) 
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L(z) ■ T(r, z), factorizing into a longitudinal profile L(z), which is a function 



of the shower depth z, and a lateral profile T(r, z), which depends on z and on the radial 
distance r from the shower center. The showers are treated as azimuthally symmetric. 



The longitudinal electromagnetic profiles are assumed to follow a Gamma distribu- 
tion, 

x <x— lg— X 

where x = fiz and z measured in units of radiation lengths Xq. a and /3 are correlated 
parameters generated using two Gaussians, whose means and widths are free parameters 
subject to tuning, and a correlation matrix hardwired in GFLASH. Longitudinal hadronic 
shower profiles are a superposition of three shower classes: 

Lhad(x) = /dep [c h Jzfh(>h) +c f J2f(xf) +aS%(xi)] , (3) 

= Xj yu^ ' Xi = $ iZi ( f = h ' f ' L ) • (4) 

Jz?h is a purely hadronic component (zt given in units of absorption lengths Ao). =£f 
accounts for the component induced by neutral pions from a first inelastic interaction 
([zf] = Xq), and J£\ originates from neutral pions occurring in later stages of the shower 
development ([z\] = Ao). Each subprofile in Eq. © is characterized by an individual cor- 
related pair (a, /3) analogously to Eq. ©. The coefficients c ; are the relative probabilities 
of the three classes expressed in terms of the fraction of showers containing a neutral 
pion (/jjo) and the fraction of showers with a neutral pion in later interactions (/J ): 

c h = 1 -fjfi , c f = fjfii} -/Jo) , a = fjfifjfi . (5) 

The global factor /d ep is the fraction of deposited energy w.r.t. the energy of the incident 
particle. When a shower is generated, correlations between all parameters are properly 
taken into account. In total, the longitudinal shower profile is described by 18 indepen- 
dent parameters for the hadronic part (the means and widths of the ct's, /3's, and of the 
fractions /d ep , f K o and f l ), and four parameters for the purely electromagnetic part. 
The lateral energy profile at a given shower depth z has the functional form 

T (r) = T^L • (6) 

The free quantity Rq is given in units of Moliere radius (for electromagnetic) or ab- 
sorption lengths (for hadronic showers), respectively. Rq is an approximate log-normal 
distribution with a mean and a variance parametrized as a function of the incident parti- 
cle energy £i nc and shower depth z- The mean is given by 

(R (E mc ,z)) = [R 1 + (R 2 -R 3 \nE inc )z} n , (7) 

where n=\(2) for the hadronic (electromagnetic) case. The spread of hadronic showers 
increases linearly with shower depth z and decreases logarithmically with E- mc . Both 
shower types have their own set of adjustable parameter values (the Ri plus three 
independent parameters for the variance), thus giving a total of 12 parameters. 

After generating the profiles, GFLASH distributes the incident particle energy in dis- 
crete interval steps following the longitudinal profile, and then deposits energy spots in 



the simulated calorimeter volume according to the lateral profile. The number of energy 
spots is smeared to account for sampling fluctuations. The visible energy is obtained by 
integrating over the energy spots and applying the relative sampling fractions Eq. O. 

4. GFLASH TUNING METHOD 

The tuning of GFLASH mostly relies on test beams of electrons and charged pions 
with energies between 8 and 230 GeV L7], and has been refined during Run II using 
the calorimeter response of single isolated tracks in the momentum range 0.5-40 GeV/c 
measured in situ. The tuning of hadronic showers follows a four-step procedure. 

Adjusting the MIP peak. Reproducing the response of minimum ionizing particles 
(MIP) is the first step since it serves as a reference for all other responses, see Eq. ([T]). 
The position and width of the simulated MIP peak is tuned to 57 GeV test beam pions 
in the electromagnetic compartments and involves the adjustment of charge collection 
efficiencies, which is handled by GEANT at this stage of the simulation. 

Setting the hadronic energy scale and shape. Next, the shape of the individual en- 
ergy distributions in the electromagnetic and hadronic compartment, the sum of both, 
and the hadronic response of MlP-like particles in the hadronic compartment are ad- 
justed. The tuning is based on 57 GeV test beam pions and involves the most sensitive 
longitudinal parameters, the means and widths of 0C\ and fi\ (steering the shower com- 
ponent induced by 7T°'s from later interactions), of /3h (which is related to the purely 
hadronic component), see Eq. ©, and of the fractions /dep an d f K o in Eqn. © and ©. 

Fixing the energy dependence. After the longitudinal profile has been adjusted at 
one energy point, the other test beam samples are employed to parametrize the energy 
dependence of the parameters, which is typically logarithmic. At this point, also the 
inclusion of in situ tracks is important in order to provide a robust extrapolation into the 
important low energy region where test beam data are not available. A functional form 
f(Ei nc ) = a + btanh(clnE mc + d) is used to describe the dependence of the fractions 
in Eqn. © and © on the incident particle energy E- mc . Distinct parametrizations in the 
central and plug calorimeter were introduced to account for their different sampling 
structure. Initial Run II tunes involved samples of single isolated tracks up to 5 GeV/c in 
the central and plug part. Recently, the energy evolution of /d ep and the relative sampling 
fractions, e/m and h/m, have been tuned using the total and MlP-like responses of in 
situ tracks up to 40GeV/c (Sec.©. 

Tuning the lateral profile. The lateral profile is treated almost independently from 
longitudinal profile details and is tuned solely with Run II data. The parameters Ri of 
Eq. 0, which have been initially adjusted to describe the profiles in the energy range 
0.5-5.0 GeV/c (using default values provided by the HI collaboration at higher energies), 
are now tuned up to particle momenta of 40 GeV/c (Sec. ED). 
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FIGURE 2. Left: Jet energy scale uncertainties due to calorimeter simulation and calibration published 
in [8]. Right: Di-jet balance j3 vs. 77 in data and simulation for transverse jet momenta 25-55 GeV/c. 



5. IMPACT ON CDF JET ENERGY SCALE 

We first report on the performance of the early Run II tuning embodied by CDF physics 
publications to date. Despite the good agreement with test beam data [7], the quality of 
the simulation used to be hard to control in the intermediate energy region, which is of 
particular importance since it constitutes a relevant fraction of the particle spectrum of a 
typical jet in CDF. For the initial tuning of the energy evolution of GFLASH parameters, 
minimum bias data providing tracks only up to 5 GeV were involved. Later checks based 
on special single track samples at higher momenta revealed a underestimation of the 
data. The level of agreement in the central was 2% for < 12 GeV, 3% for 12 to 20 GeV, 
and 4% for >20GeV BSD. Through convolution with a jet's typical particle spectrum, 
these numbers directly translate into the systematic uncertainties of the jet energy scale 
used by CDF to date (Fig.|2j left). The dominant contribution to the total uncertainties 
originate from discrepancies between the simulated and measured calorimeter response 
and from the low statistical precision of early Run II control samples. 

The inhomogeneity of the calorimeter response is accounted for by an r\ dependent 
tuning. Jet responses in the plug and wall part are re-calibrated w.r.t. the better under- 
stood central part using a correction derived from di-jet events, /3 = pj° be / pj lgger , which 
relates the transverse momentum pi of the non-central "probe" jet to the pi of the cen- 
tral "trigger" jet. A comparison of the simulated and measured di-jet balance (Fig. El 
right) shows that the tuning is reproducing many calorimeter particularities along r\ . 



6. IN SITU TUNING PROGRESS 

The situation in the central calorimeter has now substantially improved due to the 
development of dedicated single track triggers with high momentum thresholds up to 
15GeV/c. The samples collected so far (a total of over 20 M events) allow a precise 
monitoring of the simulation performance in the energy range 0.5-40 GeV. In addition, 
steadily increased single track samples selected in minimum bias events now allow a 
consistent tuning of the plug simulation up to incident particle energies of 20 GeV. 
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FIGURE 3. Simulated and measured (E/p) responses vs. track momentum in the central calorimeters 
CEM and CHA, the sum of both (TOT), and the CHA response of MlP-like particles (MIP). The vertical 
dashed lines indicate biases due to momentum dependent track analysis cuts and calculation methods. 



Measurement technique. The in situ approach employs high quality tracks which 
are well contained within the target calorimeter tower (given by the track's extrapolation 
into the electromagnetic compartment) and which are isolated within a 7x7 tower block 
around the target tower. The signal is defined as the energy deposition seen in a 2x2 
block (for PEM and CEM) or a 3x3 block (for CHA and PHA), respectively. For the 
background estimate, tower strips with the same r\ range but along the edges of a 5x5 
tower group around the target tower are used. The simulation is usually based on a 
particle gun using a controlled momentum spectrum of a mixture of pions, kaons and 
protons, plus a minimum bias generator on top for realistic background modeling. 

Absolute (E/p) response. Fig. |3] shows the comparison of the average E/p re- 
sponses between data and simulation based on the most recent tuning of the energy de- 
pendence of fdep, e/m and h/m. The basic idea is to adjust simultaneously the simulated 
total response (TOT) and the CHA response of MlP-like particles with CEM< 670 MeV 
using in situ tracks up to 40 GeV, keeping the test beam tuning at higher energies. In- 
deed, the simulated TOT and MIP (top) as well as the CEM and CHA responses (bottom) 
agree quite well with the data. The total (E / p) in the central part, which serves as an im- 
portant benchmark for the systematic jet energy scale uncertainties, is reproduced with 
1-2% precision within 0.5-40 GeV. The new tuning reflects more properly the current 
condition of the CDF calorimeter, which includes aging effects of the photomultipliers, 
and partially replaces the former conservative test beam uncertainties of 4 %. 
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FIGURE 4. Top: X 2 contours versus GFLASH parameters Ri and Q (see text) obtained from compar- 
isons of simulated and measured (E / p) profiles in CEM (left) and CHA (center) for 16-24 GeV/c particles, 
and a weighted combination of both (right). Bottom: Simulated (E/p) profiles versus data in CEM, CHA, 
and the sum of both (TOT) corresponding to the combined minimum % 2 position (highlighted point). 



(E I p) profile. An accurate adjustment of the hadronic lateral profile is important for 
two reasons. First, it controls the energy leakage out of the limited (E / p) signal regions. 
Thus, any profile mismatch between simulation and data causes a bias when tuning the 
absolute response. Second, leakage effects directly contribute to systematic uncertainties 
due to the correction of jet energies for the energy flow out of a jet cone [8]. 

For the tuning, an experimental profile is defined using the individual (E / p) responses 
of five 1x3 tower strips consecutive in 17, versus a relative r\ coordinate normalized to 
the v\ of the target tower boundaries (rj = denoting the center of the target tower). The 
availability of large single isolated track samples allows a straightforward systematic 
approach: In Eq. O, the constant Ri denotes the shower core, while Q = R2—R^ ln^nc 
fixes how the profile evolves with shower depth and incident particle energy. Since the 
electromagnetic and hadronic calorimeter compartments probe different stages of the 
average shower development, R\ and Q can be constrained using the profiles measured 
in the individual calorimeter compartments. The top of Fig.|4| shows a comparison of 
16-24 GeV profiles between data and simulation in terms of a standard x 2 estimator. 
CHA and CEM provide different contours of preferred parameter values, which helps to 
resolve the ambiguity due to the strong anti-correlation of R\ and Q by using a combi- 
nation of both (top right). Generally a good agreement between data and simulation is 
obtained (bottom). Thus, Ri has been fixed from 0.5 to 40 (20) GeV in the central (plug) 
calorimeter, and R2 and R3 can be extracted from the linear energy dependence of Q. 



7. ELECTROMAGNETIC RESPONSE 



The tuning of electromagnetic showers in GFLASH presents less difficulty and is, al- 
though important, not detailed in this report. The simulated electromagnetic scale, which 
has been set using electron test beam data, was validated in Run II at low momenta using 
electrons from J/\ff decays and at high momenta using electrons from W /Z decays . The 
simulation reproduces the measured (E/p) response with a precision of 1.7%. The un- 
certainty is dominated by a contribution of 1.6% due to electrons pointing at the cracks 
between the towers, whereas electrons well contained in the target tower account for less 
than 1% The crack response is complicated due to the presence of instrumentation 
(e.g. wavelength shifter) but can be monitored using electron pairs from Z decays. One 
electron leg is required to be well contained within a target tower and serves as an energy 
reference, the other leg is used as probe to scan the (E / p) profile along the tower up to 
the crack. The description of the crack response in GFLASH has thus been improved 
using a correction function applied to the lateral profile mapping of energy spots. 

8. CONCLUSION AND OUTLOOK 

GFLASH has been tuned to reproduce the average hadronic response in the CDF central 
calorimeter with a precision of 1-2% within the energy range 0.5-40 GeV/c. The electron 
response is reproduced with similar precision. An in situ tuning technique developed and 
re-fined in Run II has proved to be crucial to overcome past and current performance 
limits. GFLASH might also be a promising simulation tool for LHC experiments. It is 
more flexible than GEANT, it is tunable, and it demonstrates excellent CPU performance. 
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